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ABSTRACT 



Context. The near infrared (nIR)/optical counterparts of low mass X-ray binaries (LMXBs) are often observationally dim and reside 

in high source density fields which make their identification problematic; however, without such a counterpart identification we are 

unable to investigate many of the properties of LMXB systems. 

Aims. Here, in the context of a larger identification campaign, we examine the fields of four LMXB systems near the Galactic centre, 

in a bid to identify nIR/optical counterparts to the previously detected X-ray point sources. 

Methods. We obtain nIR/optical images of the fields with the ESO ~ New Technology Telescope and apply standard photometric and 

astrometric calibrations; these data are supplemented by 5pifzer-GLIMPSE catalog data. 

Results. On the basis of positional coincidence with the arcsecond accurate X-ray positions, we identify unambiguous counterpart 

candidates for XTE J 1637-498, IGRJ17379-3747, IGRJ17585-3057 and GX9+1. 

Conclusions. We propose tentative nIR counterparts of four LMXBs which require further investigation to confirm their associations 

to the X-ray sources. 

Keywords. X-rays: binaries - infrared: stars - X-rays: individuals: XTEJ1637-498, IGRJ17379-3747, IGRJ17585-3057, GX 9-1-1 



1. Introduction 

Given the main sequence or red giant nature of the companion 
stars of low mass X-ray binaries (LMXBs), they are intrinsically 
dim in the optical and near-infrared (nIR), and hence difficult to 
detect (see e.g., Charles & Coe 2003 for a review of the opti- 
cal and nIR properties of LMXBs). The intrinsic dimness of all 
these sources is further compounded by the high level of optical 
extin ction in the direction of the Galactic plane (ISch legel et al., 
11998 ). where many LMXBs are found dGrimm et ali i2002i) : 
hence observing at nIR wavelengths, where this is less pro- 
nounced, may increase the chance of a detection. Due to their 
positions in the Galactic plane or centre, where the optical/nIR 
(OIR) source density is extremely high, accurate X-ray or ra- 
dio positions are necessary to confidently identify the OIR coun- 
terparts of LMXBs and reduce the probability of chance super- 
positions. The required (sub)arcsecond positions are now more 
readily available due to missions such as Swift, XMM-Newton 
and Chandra, and when combined with the current generation 
of optical telescopes and instruments, increase the likelihood of 
counterpart discovery. 

If the LMXB is a persistent source or a transient source in 
outburst (Ivan Paradiis & McClintocla 1199 4*) emission from the 
active accretion disk may be brighter than the companion star 
and more likely to be detected at OIR wavelengths; at nIR wave- 
lengths it is even possible that a jet component will make a 
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significant contribution to the emission (e.g., iRussell & Fendeil 
2010). Broadband SED fitting of observed OIR counterparts can 
allow the various emission components (companion star and ac- 
cretion disk/jet if present) to be disentangled from each other 
and, along with spectral lines, allows the spectral type - and 
hence mass - of the companion to be constrained. When com- 
bined with radial velocity information from X-ray or optical 
light curves this allows constraints to be placed on the mass of 
the compact object with which to confirm the nature of the com- 
pact object as derived from X-ray spectral or timing properties. 
Furthermore the relationship between the OIR and X-ray flux 
can be used as a diagnostic with which to derive the compact 
object or companion star nature, and to identify the state of the 
system at the time of observations (e.g., Russell et al. 2006). 

Here we examine the fields of four LMXB systems near 
the Galactic centre, in a bid to identify nIR/optical counter- 
parts to the previously detected X-ray point sources (TablelTJ. 
ForXTEJ1637-498,IGRJ17379-3747,IGRJ17585-3057and 
GX 9-1-1 we identify unambiguous counterpart candidates on the 
basis of positional coincidence with the arcsecond accurate X- 
ray positions. We have also observed a further three LMXB sys- 
tems (1A1742-294, SLX 1744-300 and IGRJ17505-2644)* 
as part of our campaign to identify the nIR counterparts of 
high energy sources near the Galactic centre. Though we 
do not detail these obse r vation s, they confirm the results of 
IZolotukhin & RevnivtsevI (1201 lb and have been used to verify 
our reduction methods. In section |2] we introduce our observa- 
tions and reduction method while in section|3] after briefly intro- 
ducing each source, we detail the results of those observations. 
We summarise our findings in section |4] Throughout, positions 
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Table 1. X-ray and nIR positions for the LMXBs in our study 



Source 



X-ray 
RA 



Declination Err (") 



nIR 
RA 



Declination Err (") 



XTEJ1637-498' 


16 


37:02.67 


-49:51:40.6 


1.8 


16:37:02.66 


-49:51:40.0 


0.16 


IGRJI7379-37472 


17 


37:58.81 


-37:46:19.6 


1.4 


17:37:58.77 


-37:46:19.9 


0.32 


IGRJ17585-30573 


17 


58:29.85 


-30:57:01.6 


0.6 


17:58:29.83 


-30:57:01.6 


0.16 


GX9+13 


18 


01:32.15 


-20:31:46.1 


0.6 


18:01:32.15 


-20:31:46.2 


0.16 



Reference to X-ray positions: ' IStarling et alj ( l2008h . ^ iKrimm et al.l J200a) . ^ M. Revnivtsev (private communication) 



(J2000) are given with 90% confidence while all others values, 
including magnitudes, are given with Icr confidence. 

2. Observations and data analysis 

Optical (U, B, V, R, i) and nIR (J, H, K,) data were obtained with 
the ESO Faint Object Spectrograph and Camera (v.2; EFOSC2) 
and the Son of ISAAC (Sofl) infrared spectrograph and imag- 
ing camera on the 3.58m ESO - New Technology Telescope 
(NTT) over two observing runs in September/October 2008 
and March 2010 (Table|2]i. The NTT data used a dithered pat- 
tern of 3 and 9 exposures per final image in the optical and 
nIR respectively. Data were reduced using the IRAF package 
wherein crosstalk coiTection, flatfielding, sky subtraction, bias- 
subtraction and frame addition were carried out as necessary. 
The 3.9' X 3. 9' images were astrometrically calibrated against 
2MASS dStoutskie et al., 2006) or USNO-Bl.O (Monet etafl 
l2003h within the GAIA package and given positional eiTors in- 
clude a 0.16" 2MASS systematic uncertainty (TablelTJ. 

PSF photometry was carried ou t on the final images using the 
DAOPHOT package (IStetsonl[T987h within IRAF. The magnitude 
of the source of interest in each field was calculated relative to a 
number of comparison stars in the field, including the scatter as 
a measure of error. T he c omparison stars were calibrated against 
iPerssonetaP (Il998h or iLandolJ (Il992h photometric standards 
(though in the case of IGR J17379-3747, no suitable optical 
standards were observed so magnitudes were calibrated against 
USNO-Bl.O sources). If a field was observed on more than one 
night we tested for variability, using only the error associated 
with the relative magnitude. As no variability was found, the tab- 
ulated magnitudes are calculated from the weighted average rel- 
ative magnitudes and, since the exposures being compared were 
equal, the quoted exposures are of a single image on ly. The equa- 
tion,/-/ = (0.247 +0.003)(/;-/) (iJordi et al.Ll2006l) . was used to 
transform the cataloged / magnitudes of the standard stars into / 
magnitudes with which to calibrate the images. It was again used 
to transform the observed / band magnitude of the object into ap- 
propriate / magnitudes. Upper limits are approximated from the 
dimmest observable object in the region of interest. 

To supplement our own obs ervations, we utilise data from 
the Spitzer Space Telescope's dWerner et al.L l2004l) Galactic 
Legacy Infrared Mid -Plane Survey Extraordinaire (GLIMPSE; 
iBeniamin et al.l l2003h . when available. GLIMP SE was carried 
out by the IRAC instrument (iFazio et al.L l2004l) aboard Spitzer, 
spans 65° either side of the Galactic center up to ±2-4° in lati- 
tude and includes mosaic images and catalog entries at 3.6, 4.5, 
5.8 and 8.0 yum. 

The observed infrared magnitudes (Table|3]l were first con- 
verted to flux densities, Fy, at frequency v, then to flux per filter, 
F filler in units of photons cm^^s"'. This is done via Ffuter - 
1509.18896Fv (A/l/zi) where A and AA are the eff'ective wave- 
length and full width at half maximum of the filter in question. 



Table 2. Observation log 



Source 


Night of observation 


Filters 


XTEJ1637-498 


September 9 2008 


optical/nIR 




September 20, 212008 


optical, nIR 


IGRJ17379-3747 


September 29 2008 


optical/nIR 




October 2, 3 2008 


optical, nIR 


IGRJ17585-3057 


March 29 2010 


Ks 


GX9+1 


March 28 2010 


Ks 



Note that the second epoch of observations, if done, were carried 
out over 2 nights: optical on the first, nIR on the second. 



Table 3. Apparent magnitudes of the LMXB counterparts 



Source 






Filter 


Exp (s) 


Magnitude 


XTEJ1637-498t 






U 


180 


21.51 ±0.14 


B 


180 


20.77 ± 0.05 


V 


180 


18.90 ± 0.05 


R 


180 


17.66 ± 0.05 


I 


180 


16.24 ±0.16 


J 


90 


14.80 ± 0.09 


H 


90 


14.13 ±0.10 


Ks 


90 


13.55 ±0.07 


3.6fim 




12.87 ± 0.09 


4.5 ^m 




12.51 ±0.10 


5.8yum 




12.41 ±0.17 


8.0yum 




12.28 ± 0.24 


IGRJ17379-3747t 






R 


180 


21.5 ±0.3 


I 


180 


20.0 ± 0.4 


J 


90 


18.3 + 0.1 


H 


90 


> 18.0 


Ks 


90 


> 18.0 


IGRJ17585-3057 






Ks 


54 


14.18 ±0.05 


3.6/im 




11.38 ±0.05 


4.5 jum 




10.71 ±0.05 


5.8yum 




10.07 ± 0.05 


8.0yum 




9.30 + 0.03 


GX9+1 






Ks 


90 


15.35 ± 0.04 



t For observations with 2 epochs, magnitudes were consistent so we 
present the weighted average and, since the exposures being com- 
pared were equal, the exposure of a single image. 



XSPEC compatible files, for Spectral Energy Distribution (SED) 
fitting, are then produced from the flux per filter value using the 
FTOOL, flx2xsp. 
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Sg*J17585-3b5> ^^iii" • 



GX9+1 






Fig. 1. NTT images (60" x 35") of the fields of XTE J1637-498, IGR J17379-3747, IGRJ17585-3057 and GX9+1 with the X-ray 
90% position marked. IGR J 17379 -3747 is an EFOSC2 7J band image, while all others are Sofl/Ts band images; all have a common 
scaling and orientation. 



3. Results 

3.1. XTEJ1637-498 

The transient X-ray source, XTEJ1637-498, was first detected 
by the Rossi X-ray Timing Explorer (RXTE) at a 2- 1 ke V flux of 
2-4 mCrab; an outburst wh ich lasted approximately from August 
25 to September 5 2008 (iMarkwardt et al.L l2008al) . Follow up 
observations by the Swift X-ray Telescope (XRT) allowed the 
position to be refined (Table[Tl Starling et al. 2008). Based on 
a po wer law fit (photon in dex 1.5 + 0.4) of the XRT spec- 
trum, IWiinands et alj (l2008l) suggested that the source was an 
LMXB, though they cautioned that another syste m type (i.e., 
high mass X-ray binary) could not be excluded. Starling et al. 
also noted that an optical source, consistent with 2MASS object 
J16370267-4951401, was detected in the v band by the Swift 
Ultraviolet & Optical Telescope (UVOT), within the XRT posi- 
tional error. No further analysis of this source was published so 
its true classification remains unknown. 

Within the X-ray positional error, we find that an optical 
and nIR source is detected in all filters (Figure[T] best see- 
ing was ~ 0.95" in nIR filters and 1.0-1.1" in optical filters). 
From examination of our second epoch nIR images we derive 
a position of RA, Dec = 16:37:02.66 -49:51:40.0 (±0.16"), 
consistent with 2MA SS object J 16370267 -495 1401, noted by 
IStarlingetal] (120081) . We find no variabiHty (< 0.05 magnitudes) 
of the counterpart between the two epochs so present only the 
average values (TableO. We note that the observed magnitudes 
in J, H and Ks are consistent with the cataloged 2MASS values; 
this is not surprising, given that the observations were obtained 



after the X-ray outburst had faded. Furthermore we note that this 
source is also tabulated in the GLIMPSE catalog (G335.4256- 
01.7828). 

Figure|2] shows the spectral shape of the counterpart of 
XTE J1637-498, including the GLIMPSE data. Also plotted are 
the fluxes coiT ected f or a Galactic extinction in that direction 
(Schlegeleta llll998h of Eb-v = 2.14^ (A^ ~ 0.8, Ay -7.1; 
Cardelli et al.ll 19891) . though this should be considered as an up- 
per limit. Though we are unable to find a single model - for 
any value of Galactic extinction - that fits both the NTT and 
GLIMPSE data, we find that the NTT data alone can be fit by 
black body radiation at a temperature of 3100-12800K and cor- 
responding extinctions of 0.4 < Eb-v < 2.0. However, when ex- 
tended to the Spitzer data, this model underestimates the flux at 
those wavelengths. In fact, the Spitzer data itself implies a power 
law (Fy oc v") of spectral index 1 .0 < o- < 1 .7 for the same range 
of extinctions as implied by the black body fit to the NTT data. 
A power law of similar slope is not a good fit to the NTT data 
alone or the NTT and Spitzer data combined (xl » 1). This 
implies that there is excess emission at the Spitzer-IRAC wave- 
lengths (see section |4]i, either due to intrinsic IR excess or an 
excess at the time of the Spitzer observations which, as catalog 
magnitudes, were not simultaneous with the NTT observations. 



' Note that these extinctions should be treated with caution as esti- 
mates so close to the Galactic plane (< 5 deg) are unreliable. 
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Fig.2. Flux versus frequency plot for XTEJ1637-498: uncor- Fig. 3. Flux versus frequency plot for IGRJ17379-3747: uncor- 
rected (lower) and corrected (upper) for Galactic extinction in reeled (lower) and corrected (upper) for Galactic extinction in 
that direction of Eb-v - 2. 143. that direction of Eb-v = 1 -506. 



3.2. IGRJ1 7379-3747 

Outbursts from IGRJ17379-374'fl were first reported 
from the Interna tional Gamma-R a y Ast rophysics Satellite 
(INTEGRAL) by IChelovekov eTaP d2006l) with a fuller in- 
ve stigation of INT EGRAL and ot h er ob servations presented 
in IChelovekovXGrebenevl (l2010l) . Iciielovekov e t al. (2006) 
classified the source as an X-ray burster - a weakly magnetized, 
accreting neutron star in a LMXB system - on the basis of a 
type-I burst or a thermonuclear explosion, due to the unstable 
burning of H and He on the neutron star surface (iHoffman et al.L 
11978 ). During the 2008 outburst t he X-ray posit i on wa s refined 
to arcsecond precision (Tabled iKrimm et alJ l2008l) by the 
Swift-XRT 

At the X-ray position a new source is detected in our / and 
R images on both nights (Figure[Tl best seeing was ~ 1.4" in 
nIR filters and ~ 1.2" in optical filters). The source exhibits no 
variability (< 0.05, < 0.1 magnitudes respectively) between the 
2 epochs which are separated by 3 days (we present the aver- 
age magnitudes). The lack of variability is not surprising since 
observations were taken after the 20 08 September 2-21 X-ray 
outburst had faded (iMarkwardt et alil2008 b). The source is also 
detected in the second epoch J band image but not in the first, 
though the magnitude limit on that night is consistent with the 
detected magnitude. We note that the counterpart is at the edge of 
the psf of a nearby, bright 2MASS object (J17375858-3746228; 
J = 13.210,//= 12.067, /T, = 11.743). There was no detections 
in the U, B and V bands and the lack of suitable optical stan- 
dard stars on either night have lead us to omit their magnitude 
limits. The position of the optical/nIR counterpart was derived 
as RA, Dec = 17:37:58.77 -37:46: 19.9 (+0.32") from the R and 
J band images. The spectral shape of the proposed OIR coun- 
terpart of IGR J17379-3747 (Figure[3]l is not particularly con- 
straining but is consistent with stellar-like black body emission 
for a realistic range of temperatures, at e xtinctions up to the va lue 
for Galactic extinction in that direction dSchlegel et al.L Il998h of 

Eb-v = 1.50^ (A^ ~ 0.5: ICardem etanil989h . This region is 



^While IChelovekov et alj ( 120061) and |ChiiovitovXGrebini3 
1 I2OIOI) refer to the source by differing IGR designations 
(IGRJ17364-2711 and IGRJ1 7380-3749), herein we shall adhere to 
that of the INTEGRAL catalog JBird et alil2007l . l201Cl) . 
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Fig. 4. Flux versus frequency plot for IGR J 17585-3057: uncor- 
rected (lower) and corrected (upper) for Galactic extinction in 
that direction of Efi-v = 1.302. 



not covered by the GLIMPSE survey so we are unable to extend 
the SED redwards. 



3.3. IGR J1 7585-3057 



Detected by INTEGRAL (iBird et all l2007l l2010l) . 
IGRJ17585-3057 is assumed to be a LMXB based on 
Gala ctic position and the lack of an alternative classifica- 
tion dRevnivtsev et al.L 120081 Monitored as p art of both the 
INTEGRAL all-sky survey in hard X-rays dKrivonos et all 
2007^ and the 5w/ff-BAT 58-Month Hard X-ray Survey 
(Baumgartneret al., 2010) there has been no evidence of flux 
variation. No further analysis of this source was published so its 
true classification remains unknown, though an accurate X-ray 
position was obtained by Chandra (TablelT] M. Revnivtsev, 
private communication). 

At the X-ray position we find a bright, Ks - 14.183 + 0.05 
counterpart in our single SofI image of the field (Figure[T] see- 
ing ~ 0.8") and we derive a position of RA, Dec = 17:58:29.83 
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-30:57:01.6 (+0.1"). This coincides with a GLIMPSE catalog 
source (G359. 6862-03.4261), the magnitudes of which we re- 
produce in Table[3] From the SED (Figure|4|, it is clear that the 
extrapolation of the Spitzer data overestimates the Ks magni- 
tude significantly. In fact the Spitzer data implies a spectral in- 
dex of ff = -0.37 + 0.10 (when optical e xtinction is se t to tha t 
of Galactic extinction in that direction flSchlegel e t al.L Il998h . 

Eb-v = I.3O2P, A/f ~ 0.4:I Cardeni et al.ll9 89) but then drops by 
almost an order of magnitude, requiring a steepening to a > 5 to 
agree with the Ks flux. This implies either intrinsic excess emis- 
sion at the Spitzer-IRAC wavelengths or it is possible that the 
source is variable, though without multiple photometric obser- 
vations, we cannot test this. 

3.4. GX9+1 

The persistent X-ray source GX9+1, a. k.a. SgrX- 3 , was 
initially discovered by sounding rockets dBradt et all Il968t 
iMaveretal., 1970) and later confirmed by the first ded- 
icated X-ray astronomy satellite, Uhuru (as 2U 1758-20; 
iGiacconi et alj 119711 |1972|) . Thereafter, it was observed by 
most other X-ray missions which revealed that the source 
is likely an 'Atoll' source, a class of neutron star LMXB 



(e.g., Langmeier et al.l 119851 : iHasinger & van der KUs 1989: 
ISchulz et alJl989l). Previous searches for an optica l or nIR coun- 
terpart (iGottwald et all Il99ll: iNavlor et all 11991') did not yield 
any plausible candidates due to the lack of an accurate X-ray 
position, which is now available from Chandra observations 
(Table[Tl M. Revnivtsev, private communication). 

At the X-ray position we detect a. Ks - 15.35 + 0.04 source 
in our single Sofl image of the field (Figure[TJ seeing ~ 1 .4") 
and we derive a position of RA, Dec = 18:01:32.15-20:31:46.2 
(+0.16"). We note that there is also a source directly to the 
West which appears in the GLIMPSE (G009.0755-H01 . 1545) and 
2MASS (J18013196-2031464; Ks = 13.263 + 0.055) catalogs 
but we find that it is > 6cr from the X-ray position and is hence 
likely unrelated to GX 9-1- 1 . 

4. Discussion and conclusions 

To approximate the probability of a chance superposition of 
the 90% X-ray positions with random sources in the respec- 
tive fields, we calculate P ^ p^ x AsrA where p^ is the sur- 
face area number density of observed sources down to the lim- 
iting magnitude and A^^r is the area of the X-ray positional er- 
ror. The probability of chance coincidence for XTEJ1637-498 
and 1GRJ17379-3747 are quite high at 86% and 28% respec- 
tively, due to a large number of sources in the deep fields but 
more significantly, the radius of the X-ray positional errors. 
The importance of sub-arcsecond positions from e.g. Chandra 
is further illustrated by the significantly lower probabilities of 
IGRJ17585-3057 and GX9-h1 (8% and 7%) which have com- 
parable number densities to the previous two fields but much 
more accurate X-ray positions. The probability of chance su- 
perposition could be reduced further if we limited our search 
to sources of specific magnitude ranges but due to a lack of in- 
formation on these sources, we make no such assumption. 

On the basis of positional coincidence with the arcsecond 
accurate X-ray positions, we suggest unambiguous (in so far 
as there is only one per source) counterpart candidates for four 
LMXBs near the Galactic centre. We caution that given the prob- 
ability of chance superpositions in these crowded fields, the as- 
sociations are tentative and require further investigations of e.g.. 



spectral type, variability, distance, to confirm the associations 
with the X-ray sources. The spectral shapes of the counterpart 
candidates are not particularly constraining but, in general, could 
be consistent with stellar-like black body emission in the opti- 
cal/nlR but with an excess at the longer, Spitzer, wavelengths, 
as for XTE J1637-498. The excess could be due to warm cir- 
cumbinary dust, likely related t o the formation of the compact 
object (e.g.. lRahoui et alj 2010h or to non-thermal, synchrotron 
jet emission (e.g.. lGelino et al.ll2010l) . Further, possibly simulta- 
neous, photometric and spectroscopic observations are required 
to disentangle the different emission components and identify 
the spectral type of the companion star, if it is detected over the 
other possible emission mechanisms. 
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